Determination of Chiral Perturbation Theory low energy constants from a precise 
description of pion-pion scattering threshold parameters 
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We determine the values of the one and two loop low energy constants appearing in the Chiral 
Perturbation Theory calculation of pion-pion scattering. For this we use a recent and precise sum 
rule determination of scattering lengths and slopes that appear in the effective range expansion. 
In addition we provide new sum rules and the values for these coefficients up to third order in the 
expansion. Our results when using only the scattering lengths and slopes of the S, P, D and F waves 
are consistent with previous determinations, but seem to require higher order contributions if they 
are to accommodate the third order coefficients of the effective range expansion. 

PACS numbers: 12.39.Fe,13.75.Lb,11.55.Hx 



I. INTRODUCTION 

The smallness of the u and d quark masses together 
with the spontaneous SU(2) chiral symmetry breaking 
of QCD, which implies the existence of the correspond- 
ing Goldstone Bosons — the pions — allow us to write a 
low energy effective theory for QCD, organized as a sys- 
tematic expansion in pion masses and momenta. This 
is known as Chiral Perturbation Theory (ChPT) [J, 0, 
built as the most general low energy expansion of a La- 
grangian containing just pions, which is compatible with 
the symmetry constraints of QCD. In particular, the first 
order is determined by the scale of the spontaneous sym- 
metry breaking, identified with the pion decay constant 
/tt, and the pion mass M n . The expansion is then carried 
out in powers of p 2 /(47r/ 2 ), where p denotes generically 
either the pion mass or momenta. The details of the QCD 
underlying dynamics are encoded in a set of low energy 
constants (LECs) , which multiply the independent terms 
that appear in the Lagrangian at higher orders. Note 
that all loop divergences appearing in a calculation up 
to a given order can be reabsorbed by renormalization of 
the LECs up to that order. In this process the LECs ac- 
quire a dependence on the renormalization scale fj,, which 
is canceled with that present in the loops. In this way 
calculations are rendered finite and scale independent to 
any given order of the expansion. 

Only certain combinations of LECs appear in ttt: scat- 
tering up to a given order. As we commented above, to 
leading order 0(p 2 ) there are no LECs. Within the SU(2) 
formalism, to next to leading order (NLO), or 0(p 4 ), 
which corresponds to a one-loop calculation, only four 
LECs, called li,h,l3 and U appear in the amplitude, al- 
though two of them ^3 and Z4, do so through the quark 
mass dependence of the pion mass M T and decay con- 
stant /tt. To next to next to leading order (NNLO), or 
0(p 6 ), six possible independent terms appear |3| multi- 
plied by six constants, hi with i — 1...6, which can be 
reexpanded in powers of the pion mass in terms of the 
four one-loop Ik and six new NNLO LECs, denoted by 



n 0. After renormalization these constants develop a 
scale dependence becoming l\ (/x) and r£(/i), whereas the 
bi remain renormalization scale independent. 

Concerning the 0(p A ) LECs, we refer the reader to 
Q for a recent compilation of lattice QCD and to Q for 
some other estimates from quark-model-like calculations. 
It is also worth noticing that the bulk of their values can 
be explained by the effect of integrating out heavier res- 
onances and, actually, seems to be saturated by the vec- 
tor multiplets Q, plus a small contribution from scalars 
above 1 GeV, and an additional sizable contribution from 
kaons in the case of the SU(2) formalism [8]. Some esti- 
mates from resonance saturation have also been obtained 
for the 0(p 6 ) parameters 0]. However, since perturba- 
tive QCD cannot be applied at very low energies, it is 
particularly difficult to obtain the values of these LECs 
from first principles and, with few exceptions, the LECs 
have been determined best from the comparison with ex- 
periment 

Our aim in this work is to use a very recent dispersive 
analysis of data in 13(, which includes, among others, 



the latest very precise and reliable results on K e ± decays 
from the NA48/2 collaboration^ . in order to determine 
the values of the 0(p A ) and 0(p 6 ) LECs that appear m 
the 7T7T scattering amplitude. Since ChPT is a low en- 
ergy amplitude, obtained as a truncated series in powers 
of (p I '4-7T j V) 2 ™ , we will compare with data at threshold. 
In particular, we will obtain the LECs from fits to the co- 
efficients of the momentum expansion of the amplitude 
around threshold, usually known as the effective range 
expansion. The coefficients of this expansion, even up to 
third order, are also becoming reachable for lattice calcu- 
lations, although still only limited to the highest isospin 
channels 

Thus, in the next section, after introducing the nec- 
essary notation, we present the experimental determi- 
nation made in [l3[ of the threshold parameters up to 
second order of the effective range expansion. The most 
precise way to obtain these parameters is by means of 
the Froissart Gribov representation and other sum rules, 
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which we briefly review. Moreover, in this work, we de- 
rive a new set of sum rules to calculate the third order 
threshold parameters, including their uncertainties, up 
to the F waves. Some of these third order parameters 
are of relevance to obtain the values of the LECs, since 
their leading order contribution is directly proportional 
to combinations of LECs. Actually, we carefully explain, 
for each threshold parameter, what is its leading order 
ChPT, and from what part of the calculation stems from. 
In Sect. IHII we first perform a fit of some of these param- 
eters using just the 0(p 4 ) ChPT result, paying particular 
attention to an estimate of the systematic uncertainties 
in the parameters, which is of relevance for the role they 
will play in the determination of different LECs. Still, 
just by trying to describe a few threshold parameters, 
we are able to show that the 0{p A ) approximation is not 
enough to describe the data at the present level of preci- 
sion. In Sect. IIVI we determine the best hi constants and 
LECs that appear in the two-loop calculation. We will 
show that one can obtain only a relatively fair descrip- 
tion in terms of bi parameters, although the fact that the 
\ 2 /d.o.f. of the fits is somewhat larger than one, sug- 
gests that at the present level of precision, higher order 
contributions seem to be required. In Sect.fVl we briefly 
discuss and summarize our results. 



II. THRESHOLD PARAMETERS 
A. Notation 
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TABLE I: Contribution to threshold parameters from differ- 
ent orders and kinds of terms within ChPT, as explained in 
the text. Recall that, due to Bose symmetry, those for S and 
D waves may have either isospin or 2, whereas those for P 
and F necessarily have isospin 1. 

them, generically, threshold parameters. Let us remark 
that it is usual to provide the values of these parameters 
in units of M w , and so we will do in what follows. In 
addition, for odd waves we will drop the isospin index, 
since it can only be / = 1 due to Bose symmetry. Finally, 
we will also make use of the standard spectroscopic nota- 
tion, where the I = 0, 1, 2, 3... are called S,P,D,F... waves, 
followed by the value of the isospin. As it is customary, 
when the angular momentum is odd we omit the isospin, 
since it is fixed to one by Bose symmetry. 



The amplitude for 7T7t scattering is customarily decom- 
posed in terms of partial waves t\ % of definite isospin / 
and angular momentum £, as follows: 

F 7 (M) = -V(2*+l)tf(«)P<(cos0), (1) 

7T ' — ' 

I 

1 f 1 

4(s) = — J T I (s,t,u)Pe(cos0)d(co86), (2) 

9 being the scattering angle, Pi the Legendre polyno- 
mials, s, t, u the usual Mandelstam variables satisfying 
s + t + u = 4:M%, and T stands for the amplitude. In the 
elastic regime, the partial waves are uniquely determined 
by the phase shifts Sj as follows: 



(3) 



where a = 2p/^/s = y/l — AM%/s and p is the CM mo- 
mentum. With this normalization, the effective range 
expansion of the real part of a partial wave can be writ- 
ten as: 

-^-Re4(s) =p^(a n + b u p 2 + ^c elP i + ...), (4) 

where the an are usually called scattering lengths, the 
ba slope parameters, the ca shape parameters, and all of 



B. Structure of ChPT calculations 

At this point it is relevant to discuss how the different 
orders of ChPT contribute to each threshold parameter 
studied here, which we have gathered in the first column 
of TablelH Let us start with the leading order, 0(p 2 ), of 
the ChPT amplitude, which is a first order polynomial 
in terms of Mandelstam variables s, t, u and M 2 . Since s 
is independent of 9 whereas t and u are first order poly- 
nomials in cos (9, the LO ChPT can only contribute to 
the a coefficients of the S and P waves and the b coef- 
ficients of the S waves, but nothing to any other wave. 
This corresponds to the second column in Table [J] 

If we now consider the 0(p A ) amplitude, we find two 
kinds of terms. First, a polynomial, which includes the 
li(fi) LECs and contains up to two powers of cos 9, so that 
it contributes to the a coefficients of S, P and D waves, 
the b coefficients of the S and P waves as well as the c 
coefficients of the S waves. This is the third column in 
Table |H However, there is another kind of 0(p 4 ) contri- 
butions, which we have represented in the fourth column. 
These come from the loop functions, called J(q 2 ), with 
two intermediate pions exchanged in any channel. These 
loop functions carry a non-polynomial dependence on t 
and u and therefore contribute to all waves, but note that 
they do not depend on ^(/i). This is the fourth column. 
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labeled "J", in Table H 

Next, to two loops, 0(p 6 ), we find three kinds of terms: 
First, pure polynomial terms containing the LECs, 
which can contribute to the a coefficients up to the F 
wave, b coefficients up to D waves, and c coefficients up 
to P waves. These appear in the fifth column of Table HI 
In addition there are terms contributing to all waves, as 
shown in column six, containing a single one-loop func- 
tion and k(n) LECs. Finally, there are also terms without 
LECs, which correspond to the last column, containing 
either two one-loop functions or one two-loop function, 
that we generically call K(q 2 ). 

Therefore, we see that only the as,bs and ap have 
a leading order contribution independent of the LECs. 
A priori, one could expect that the best observables to 
determine the li are those whose leading term is 0(p 4 ), 
as it is the case of the aoo and ap>2, which, actually, 
have been frequently used to determine the value of l± 
and I2 [2|. Nevertheless, let us remark that there are still 
other threshold parameters whose leading contributions 
are proportional to a combination of U, namely, bp, cso 
and Cs2- However, these are much harder to determine 
reliably from experiment and had not been used so far. 
For this reason, in the next subsection we will provide 
two new sum rules for their determination. To extract 
the 0(p 6 ) contributions and the r$ LECs is more difficult, 
not only because they are generically an smaller effect, 
but also because they appear together with 0(p 4 ) terms 
or with 0(p 6 ) terms containing the U and a one-loop 
function. 



C. Threshold parameters from sum rules 

As commented above, we are going to use the recent 
and simple data parametrizations obtained in [131 ], The 



relevance of those parametrizations is that they are ob- 
tained from data fits which have been highly constrained 
to satisfy three sets of dispersion relations within uncer- 
tainties: Forward dispersion relations up to 1420 MeV, 
and Roy equations as well as once subtracted Roy-like 
equations up to 1100 MeV. Above 1420 MeV, Regge ex- 
pressions, assuming factorization, where fitted to NN , 
Nit and 7T7t total cross sections, and used inside the inte- 
grals, allowing for the variation of the Regge parameters 
within the constrained fits to data. In that work, the 
values of the a and b threshold parameters were also pro- 
vided for the SO, S2, P, DO, D2 and F waves, namely, all 
the combinations of / = 0,1,2 and £ — 0,1,2,3 allowed 
by Bose symmetry when considering pions as identical 
particles. With the aim of minimizing the uncertainties, 
they were obtained from sum rules, with the only excep- 
tion of the 5aso + 2as2 combination, which is orthogo- 
nal to the one appearing in the Olsson sum rule. The 
results from sum rules were consistent with, but more 
accurate than, those directly obtained from the simple 
phenomenological parametrizations. Let us nevertheless 
remark that the sum rules used in [l3[ as well as those 
we will describe below have a very small dependence on 
the high energy region. 

These results provided us with 12 observables deter- 
mined from experiment, which we list in Table HQ that 
we want to fit using four Zj LECs in the one-loop case 
and six parameters bi in the two-loop case, which can 
be parametrized in terms of ten LECs. Moreover, in or- 
der to enlarge the set of observables that we include in 
our fit, we will also provide here the calculation of the 
third order coefficient c of the effective range expansion 
in Eq.Q above. These are five additional observables. 

For this purpose we will use the Froissart-Gribov sum 
rules, which, for I > allow us to write the c parameters 
as: 
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where the primes denote the derivative with respect to 
cos 6. These formulas allow us to calculate the c parame- 
ters for the P, DO, D2 and F wave, that we list in Table ILT1 
Note that the resulting values are all very accurate with 



the exception of the cp coefficient, and that the above 
sum rules are not applicable to the scalar case. These 
are the reasons why we provide here three new sum rules, 
one for cp and two for csq and cs2. 
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Let us note that, in all these sum rules, we have written 
several terms together inside square brackets to empha- 
size that they do not converge separately. The derivation 
is similar to the sum rules obtained for bp, bso and &S2 in 
[IH [l6j . They correspond to the threshold limit, taken 
from above, of the second derivative of a forward dis- 
persion relation for the F I=1 , F 0+ and F 00 amplitudes, 
respectively. Let us recall that F 0+ = F I=2 /2 + F I=1 /2 
whereas F 00 = 2F I=2 /3 + F I=0 /3. 

In the above sum rules we have explicitly got rid of the 
principal part that appears in the dispersion relation by 
using: 

r°° dx 

PP 7„ F^vS^^ ' (9) 



III. 0(p 4 ) FITS 

Before presenting the fits to the full two-loop ChPT 
results, it is instructive to try to fit the threshold param- 
eters by means of the one-loop ChPT amplitudes. This 
will help us check the stability of the LECs values and 
the need for higher order counterterms, but it will also 
help us illustrate our different fitting strategies in order 
to deal with systematic uncertainties. 

Let us recall once more that to 0(p 4 ) only four LECs, 
appear in tttt scattering, customarily denoted by Zi, I4,, 
which are basically the at the p, = M v scale and nor- 
malized so that they have values of order one 0. Note, 
however, that I3 and Z4 only appear through the quark 
mass dependence of M„ and f T , respectively, and there- 
fore we cannot expect much sensitivity to these two pa- 
rameters from fits to the coefficients of the momentum 
expansion of amplitudes. 

In addition, from Table U we see that, up to 0(p 4 ), 
only ten observables carry any dependence on the LECs: 
for half of them, aso, as2, ip, &so an d 652, the lead- 
ing contribution is 0(p 2 ), whereas for the other five the 
leading contribution is directly of 0(p 4 ). Therefore, we 
expect the latter to be more sensitive to the LECs, but 
also to the higher order corrections that we are neglect- 
ing. 





CFD 


Sum rules 


Best value 




0.221 ± 0.009 




0.220 ± 0.008 


a S2 (xl0 2 ) 


-4.3 ±0.8 




-4.2 ±0.4 


2aso — 5as2 


0.657 ± 0.043 


0.648 ±0.016 


0.650 ±0.015 


a P (xl0 3 ) 


38.5 ± 1.2 


37.7 ± 1.3 


38.1 ±0.9 


aoo(xl0 4 ) 


18.8 ± 0.4 


17.8 ±0.3 


17.8 ±0.3 


a D2 (xl0 4 ) 


2.8 ± 1.0 


1.85 ±0.18 


1.85 ±0.18 


of(x10 5 ) 


5.1 ±1.3 


5.65 ±0.23 


5.65 ±0.23 


bso 


0.278 ± 0.007 


0.278 ± 0.008 


0.278 ± 0.005 


MxlO 2 ) 


-8.0 ±0.9 


-8.2 ±0.4 


-8.2 ± 0.4 


6p(xl0 3 ) 


5.07 ±0.26 


6.0 ±0.9,5.48 ±0.17 


5.37 ±0.14 


&do(x10 4 ) 


-4.2 ±0.3 


-3.5 ±0.2 


-3.5 ±0.2 


&D2(X10 4 ) 


-2.8 ±0.8 


-3.3 ±0.1 


-3.3 ±0.1 


6f(x10 5 ) 


-4.6 ±2.5 


-4.06 ±0.27 


-4.06 ±0.27 


cso(xl0 2 ) 


-0.12 ± 1.22 


0.7 ±0.8 


0.45 ± 0.67 


c S2 (xl0 2 ) 


3.6 ±1.8 


2.79 ±0.24 


2.80 ±0.24 


c P (xlO 3 ) 


1.41 ±0.19 


2.3 ±0.8,1.35 ±0.15 


1.39 ± 0.12 


C D0 (X10 4 ) 


5.6 ±0.4 


4.4 ±0.3 


4.4 ±0.3 


c D2 (xl0 4 ) 


5.5 ±1.6 


3.6 ±0.2 


3.6 ±0.2 


cf( xlO 5 ) 


11 ±9 


6.9 ±0.4 


6.9 ±0.4 



TABLE II: Values of threshold parameters obtained in Ref. 
[l3| together with those obtained here for the c parameters. 
The "CFD" column lists the values as obtained directly from 
the Constrained Data Fits provided in [l3|| . We also provide 
the values obtained from sum rules. We typically consider 
these our best results, except in cases when the CFD are 
competitive and not very correlated with the sum rule. Note 
that for bp and cp waves we provide two values. For cp, the 
first one, less accurate, corresponds to the Froissart-Gribov 
sum rule in Eq.([S} and the second one to the sum rule in 
Eq.©. Similarly, for bp, the first, less precise result, is from 
the Froissart-Gribov sum rule, and the second from a fast 
convergent sum rule, as explained in [l5l. [Tr^]. 

Thus, in Table ITU we show the results of our fits. First, 
we have fitted only the observables which have a lead- 
ing 0(p 2 ) contribution, since, in principle these might 
be more stable under the higher order corrections. The 
fit comes out with relatively low \ 2 /d.o.f.. Next we have 
presented a determination of l\ and l 2 , which are, in prin- 
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ciple, fixed from aoo and ao2 alone, which are not in- 
cluded in the previous fit. It is evident that the resulting 
values from those two fits are incompatible, particularly 
l\. The incompatibility is even worse when fitting simul- 
taneously the ten observables that depend on k to 0(p A ). 
These results imply that, as is well known, to the present 
level of precision the one-loop ChPT formalism is not 
enough and calls for higher order corrections. 

For instance, the effect of higher order corrections can 
be seen by fitting to the one-loop amplitude but replacing 
fir by fo in the 0(p A ) terms, since the resulting expres- 
sion is also correct up to 0(p 4 ), only differing in higher 
order contributions. This we show in row 5 of Table 11111 
Surprisingly the x 2 fd.o.f comes somewhat lower, but the 
values of the LECs come out rather different from the 
previous calculation. 

Of course, one could always try to include the effect of 
higher orders into a systematic uncertainty of the LECs, 
if one still wants to use the relatively simple 0(p ) ap- 
proximation instead of the full two-loop amplitude, at 
the expense of accuracy. In such case we propose to take 
the weighted average of the two previous fits, including a 
systematic uncertainty to cover the LECs values of both 
fits 1 . In Table lTVl we compare the resulting threshold pa- 
rameters obtained using this averaged set with the Best 
results of Table HH which we repeat under the "Data anal- 
ysis" column. We can see there that, thanks to the larger 
uncertainty, the threshold parameters obtained are com- 
patible within errors with the experimental values, except 
for bso and bp, which differ by more than three and two 
standard deviations respectively. 



IV. 0(p 6 ) FITS 

As already commented in the introduction, the thresh- 
old parameters can be described in ChPT at 0(p 6 ) 
in terms of six low energy constants, usually denoted 
,&6- Let us remark, however, that the first four 
can be separated in two parts with different chiral order, 
namely, k = b^+Ab,, i = 1,2, 3, 4, where b[ 0) = 0(m°) 

and Abi — O(rn^). The bf ^ parameters contain combi- 
nations of the four 0(p 4 ) LECs 1$, but not of the 0(p 6 ) 
LECs. In contrast, six linear combinations of the latter 
appear inside the Abi for * = 1...4 as well as in 65 and 
be, and are accordingly denoted by rj, with i = 1...6. 
Due to this 0(rn^) part in the parameters, the calcu- 
lations using the bi have an extra 0(p 8 ) piece which is 
not present when using U and Ti (or making the separa- 
tion bi = of + Abi explicit). Of course, since this is a 
higher order contribution, both descriptions are formally 



We have weighted each fit by the square of its \ 2 /d.o.f., whereas 
the uncertainty is obtained by adding in quadrature the averaged 
statistical uncertainty to the maximum difference between the 
resulting central value and the central value of each fit. 



equivalent up to 0(p 6 ). Nevertheless, there could be rel- 
evant numerical differences and, what is more important 
to us, in one case one should determine only 6 parame- 
ters, whereas in the other case there are 10 parameters 

Thus, when using k and 7"j, we may obtain spurious 
solutions or, in general, less stable values than when using 
just the six bi. That is why in this section we have decided 
to use the bi set in our fits. For completeness we provide 
in the appendix a detailed account of our results when 
parametrizing the ChPT series in terms of U and r,-. 

Therefore, and similarly to the 0(p 4 ) case, we first fit 
the ten threshold parameters as, ap, an, bs, bp, and cs, 
which have a non-zero 0(p 4 ) polynomial contribution, 
since we expect these to be more stable under higher or- 
der ChPT corrections. In the first row of Table [V] we 
show the resulting bi for this fit, which describes fairly 
well the fitted observables with a x 2 /d.o.f . — 1.2. How- 
ever, the threshold parameters which are not fitted, are 
not so well described with this set of LECs. 

Actually, when fitting all 18 threshold parameters, 
we obtain somewhat different LECs, which are shown 
in the second row of Table |V] Although not dramati- 
cally incompatible with those of the first row, we see 
differences around two standard deviations for 63 and 
64 and around three standard deviations for b\ and 65. 
Unfortunately, this second fit comes out with a rather 
poor y 2 /d.o.f . — 5.2. Therefore, it seems that we can- 
not describe all observables simultaneously with two-loop 
ChPT within the present level of precision. Higher order 
contributions seem to be required. 

Nevertheless, we have noticed that cp alone con- 
tributes almost to one third of the total x 2 ■ This might 
indicate that cp receives important higher order contri- 
butions that are not being taken into account in the 
0(p 6 ) calculation. Once again we can obtain a crude 
estimate of the size of higher order ChPT corrections, 
by changing / ff by fo in the last term of the ChPT ex- 
pansion. By far it is cp the one that suffers the largest 
change, by almost 80%. Certainly it looks as a good 
candidate to receive very large higher order ChPT cor- 
rections. 

Thus, we proceed to fit again all threshold parameters 
except cp, and the result is shown in the third row of 
Table [V] The fit quality improves sizably, but we still 
get a high y 2 /d.o.f . = 2.9, which indicates that the two- 
loop calculation may not be enough to describe even the 
remaining threshold parameters with their current level 
of precision. 

As a final attempt, we can see the effect of higher or- 
der corrections by making a fit replacing / w by fo in the 
0(p 6 ) terms, since the resulting expression is also correct 
up to 0(p e ). We show the results (without including cp) 
in the fourth row of Table fVl Surprisingly, we now obtain 
a good x 2 /d.o.f . = 1.0 , and all LECs are less than two 
standard deviations from those obtained by fitting only 
the threshold parameters which have an 0(p i ) polyno- 
mial part. We therefore conclude that, by excluding cp, 
the two-loop fit still shows some tension but by conve- 
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Fit to 


h 


h 


^3 h 


X 2 /d.o.f. 


as, bs, ap 


l.lil.O 


5.1±0.7 


-1±8 7.1±0.7 


0.23 


'ID 


-1.75±0.22 


5.91±0.10 


— — 





cs 


-2.4±0.9 


4.8±0.4 







as, bs, a P , a D , cs, bp 


-2.06±0.14 


5.97±0.07 


-5±8 7.1±0.6 


7.9 


as, bs, a P , a D , Cs, b P , using 


fo -1.06±0.11 


4.6±0.9 


0±6 5.0±0.3 


7.06 


Estimate 0(p 4 ) 


-1.5±0.5 


5.3±0.7 


-3±7 6. Oil. 2 





TABLE III: 0(p 4 ) fits to different sets of threshold parameters containing polynomial 0(j> 4 ) contributions. Note that the results 
of the three first lines are rather incompatible with each other. This is also illustrated by the large \ 2 /d.o.f. when fitting all 
the observables simultaneously. We also show two versions of such a fit, either using f^ or fo in the last order or the expansion. 
Finally, we provide an estimate of how much one should enlarge the uncertainties of the LECs if, for simplicity, one still insists 
in using the one-loop formalism. Beyond that accuracy a two loop formalism is called for. 



niently using fo the last term of the ChPT expansion, it 
can give an acceptable description of the rest of threshold 
parameters. 

For this reason, we have once more made a weighted 
average of the two fits (the one using /„• and the one 
using fo) adding systematic uncertainties to cover both 
sets. This we show in the fifth row of Table [V] where we 
can see that they are also quite compatible with previous 
determinations in the literature [Tlj | . Let us remark that, 
as emphasized in (Tlj . "the error bars only indicate the 
noise" seen in their evaluation. This would correspond 





Estimate 


Estimate 


Data 




0( P 4 ) 


o( P 6 ) 


Analysis 


aso 


0.214±0.009 


0.230 ±0.014 


0.220 ± 0.008 


a S2 (xl0 2 ) 


-4.4±0.3 


-4.3±0.4 


-4.2 ±0.4 


a P (xl0 3 ) 


38.7±1.2 


39.0±0.8 


38.1 ±0.9 


a D o(xl0 4 ) 


15±3 


16.9±0.9 


17.8 ±0.3 


a D2 (xl0 4 ) 


1.3±1.0 


1.7±0.3 


1.85 ±0.18 


a F (xl0 5 ) 




4.6±0.5 


5.65 ±0.23 


bso 


0.255±0.011 


0.271±0.007 


0.278 ±0.005 


6 S2 (xl0 2 ) 


-8.2±0.5 


-8.4±0.2 


-8.2 ± 0.4 


6p(x10 3 ) 


4.4±0.5 


5.2±0.2 


5.37 ±0.14 


&do(x10 4 ) 




-3.6±0.8 


-3.5 ±0.2 


b D2 {xW i ) 




-3.1±0.4 


-3.3 ±0.1 


&f(x10 5 ) 




-3.4±0.3 


-4.06 ±0.27 


cso(x10 2 ) 


2.3±1.4 


1.3±0.6 


0.45 ± 0.67 


c S2 (xl0 2 ) 


3.4±0.7 


2.78±0.16 


2.80 ±0.24 


c P (xl0 3 ) 




0.3±0.2 


1.39 ±0.12 


cdo(x10 4 ) 




3.6±0.2 


4.4 ±0.3 


c C2 (xl0 4 ) 




3.2±0.2 


3.6 ±0.2 


c F (xl0 5 ) 




5.4±0.4 


6.9 ±0.4 



TABLE IV: Values of the threshold parameters obtained from 
0(p 4 ) ChPT and 0(p 6 ) ChPT using the averaged sets of LECs 
from the sixth row of Table IIIII and the fifth row of Table IVl 
We also show in the third column the best values obtained 
from the data analysis given in Table [TT1 



to our uncertainties in the "All but cp, using f n row, 
whereas the error bars we provide for our final result also 
contain some crude estimate of higher order uncertain- 
ties. 

Finally, the values obtained for the threshold param- 
eters using this averaged set of LECs are shown in the 
second column of Table ITVl where we can see that, with 
the exception of cp , they are rather compatible with the 
experimental determination. 



V. SUMMARY AND DISCUSSION 

In this work, we have determined the low energy con- 
stants of SU{2) Chiral Perturbation Theory (ChPT) at 
one and two loops from the determinations of the thresh- 
old parameters obtained from sum rules using a recent 
and precise dispersive analysis of data [lH, together six 
additional observables that we have studied here. 

Threshold parameters are denned as the coefficients 
of the effective range expansion of 7T7T scattering partial 
waves, which in this work we have studied up to angular 
momentum £ = 3, i.e., S,P, D and F waves, and all possi- 
ble isospin states I = 0,1,2. The coefficients of the two 
first orders, namely the scattering lengths an and slope 
parameters bti, were already obtained from a dispersive 
analysis of data in (l3| . In addition, we have provided 
here three new sum rules to estimate the third order co- 
efficients C(j, thus adding six new observables to form a 
total set of eighteen. We have briefly reviewed how the 
different terms and low energy constants of ChPT con- 
tribute to each one of these threshold parameters. 

We have then proceeded to fit these observables, first 
within one- loop ChPT 0(p 4 ), and then to the full two- 
loop 0(p 6 ) calculation [4j. We have checked that the one- 
loop formalism is clearly insufficient to accommodate the 
present level of precision. There is a clear improvement, 
in terms of x 2 /d.o.f . when using the two-loop expansion, 
although it is still not sufficient to get a good quality fit. 
This suggests that even higher order ChPT contributions 
may still be required to describe all these observables 



7 



Fit to 


5i 


b 2 


h 


b A 


h 


b fi 


X 2 /d.o.f. 




as,b s ,ap,ao,c s ,bp 


-14±4 


14.6±1.2 


-0.29±0.05 


0.76±0.02 


0.1+1.1 


2.2±0.2 


6.0/(10-6+1)= 


=1.2 


All 


-2+3 


14.2±1.0 


-0.39±0.04 


0.746±0.013 


3.1±0.3 


2.58±0.12 


67/(18-6+1)= 


=5.2 


All but cp 


-6±3 


15.9±1.0 


-0.36±0.04 


0.753±0.013 


2.2±0.4 


2.44±0.12 


34.9/(17-6+1): 


=2.9 


All but cp, using fo 


-12±3 


13.9±0.9 


-0.30±0.04 0.726±0.013 1.0±0.3 1.93±0.08 12.5/(17-6+1)= 


=1.04 


Estimate 0(p 6 ) 


-10.5±5.1 


14.5±1.8 


-0.31±0.06 


0.73±0.02 


1.3±1.0 


2.1±0.4 






Ref. [Ill 


-12.4±1.6 


11.8±0.6 


-0.33±0.07 


0.74±0.01 


3.6±0.4 


2.35±0.02 







TABLE V: 0(p 6 ) fits to different sets of threshold parameters. In the first row we only fit to observables containing polynomial 
0(p 4 ) contributions. Note the improvement of the 0(p 6 ) description versus the 0(p 4 ) one by comparing the \ 2 /d-o.f. here 
with the corresponding one in Table IIIII Next we show the results of the fit to all the threshold parameters obtained in this 
work. Note that the fit quality is rather poor. However, most of the disagreement is caused by a single observable cp. When 
this observable is omitted, the resulting fits are of much better quality, particularly good when using fo instead of f n in the 
last term of the ChPT expansion. We also provide an estimate of the LECs uncertainties from the fits to all observables except 
cp, as a weighted average of the fits using fo or f^. Within our uncertainties, the resulting values of the bi parameters are 
very consistent with previous determinations, listed in the last row. Let us remark that, as emphasized in [ll[, "the error bars 
only indicate the noise" seen in their evaluation. This would correspond to our uncertainties in the "All but cp, using fo" row, 
whereas the error bars we provide for our final result also contain some crude estimate of higher order uncertainties. 



simultaneously. 

However, we have been able to identify that the largest 
incompatibility is due to the cp parameter. This may not 
come as a big surprise, since the largest contribution to 
the value of the sum rule that determines cp is given 
by the p resonance and its sharp rise before 770 MeV, 
which cannot be reproduced by the perturbative ChPT 
series. Actually, we have estimated, by changing /„• from 
its physical value to its value in the chiral limit in the 
last term of the 0(p 6 ) expansion, that this observable is 
a natural candidate to receive very large corrections form 
higher ChPT orders. 

Hence, if cp is omitted, the quality of the two loop 
fit improves, although there is still some tension in the 
parameters to describe the remaining threshold parame- 
ters. Nevertheless, by conveniently using /o in the last 
term of the ChPT expansion, the two -loop expansion can 
provide an acceptable description of the rest of threshold 
observables. The ChPT parameters thus obtained, for 
which we provide statistical as well as an estimate of sys- 
tematic uncertainties, are fairly compatible with previous 
determinations. 
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Appendix A: Fits to U and fi 

As commented in Section UVT the two- loop 7T7t scatter- 
ing amplitudes can be recast in terms of six independent 
terms multiplied by their corresponding low energy con- 
stants bi. In turn, these bi can be rewritten in terms of 
the four 0(p 4 ) LECs that appear in the Lagrangian and 
six combinations r; of the 0(p 6 ) LECs. The difference 
between writing the amplitude in one way or the other 
is 0(p s ). However, despite increasing the number of pa- 
rameters to ten, the 0(p e ) amplitude still provides just 
six independent structures. As a consequence, the fits 
in terms of U and are much more unstable, and can 
even lead to spurious solutions. For this reason we have 
explained the fits in terms of 6, in the main text, and we 
have relegated the U, ri fits to this appendix. 

Let us then revisit the fits of Table |V] where we fit all 
the observables, or all but cp, but recasting the ampli- 
tudes in terms of U and f^. 2 The resulting values are given 
in Table IVT1 We observe the same pattern as before: the 
fit to all parameters still has a rather high \ 2 /d.o.f. = 5, 
because, although the total % 2 has decreased from 67 to 
42, the number of degrees of freedom has increased by 4. 
Let us remark that f\ and r% have central values many 
orders of magnitude bigger than expected, but their un- 
certainties are comparably large. This means that we do 
not have any real sensitivity to these parameters. 

Once again, the largest contribution to the x 2 is due 
to cp, and thus we remove it from the fits, as we did in 
the main text. When so doing, the ^/d.o.f. becomes 
less than one, yielding a statistically acceptable fit, but 
of course, the uncertainties are still huge for f\ and fi 
and very large for ^3, and r^. The central value of Z3 



2 We give the fit results in terms of fi, defined as fi = (47r) 4 r;, 
which have a size of order one. 
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Fit to: 


h 


h 


h 


h 


fi 


f 2 


^3 


fi 




fa 


X 2 /d.o 


/■ 


All 


-0.88±1.43 


5.1± 0.8 


-49±10 


4.5±1.3 


-984±335 


-101±302 


-5.7±26 


-13±1S 


1.6±0.9 


0.45±0.33 


42 


= 5 


18-10+1 


All but cp 


-2.2±1.5 


5.6±0.8 


-20±11 


10±2 


276±845 


-1361±549 


34±36 


-38±19 


0.67±0.94 


0.66±0.34 


6.7 


= 0.8 


17-10+1 


All but cp with /o 


-0.5±1.0 


4.2±0.6 


-6±8 


6.6±1.1 


46±450 


-356±238 


4±15 


-9±9 


1.5±0.6 


0.5±0.2 


4.7 


= 0.6 


17-10+1 


Constrained fit to: 


All but cp 


-0.11±0.16 


4.2±0.1 


3.3 


5.8±0.4 


-1.5 


3.2 


-4.2 


-2.5 


3.1±0.5 


0.85±0.15 


68 
17-5+1 


= 5.2 


All but cp with /o 


0.5±0.2 


3.9±0.1 


3.3 


5.1±0.3 


-1.5 


3.2 


-4.2 


-2.5 


1.4±0.4 


0.47±0.12 


15.7 
17-5+1 


= 1.2 


Estimate 0(p 6 ) 


0.4±0.5 


3.9±0.3 


3.3 


5.2±0.7 


-1.5 


3.2 


-4.2 


-2.5 


1.7±1.5 


0.5±0.3 







Estimate 0(p 4 ) -1.5±0.5 5.3±0.7 -3±7 6.0±1.2 

Refs.[4, jj, Ft] -0.4±0.6 4.3±0.1 3.3±0.7 4.4±0.2 -1.5 3.2 -4.2 -2.5 3.8±1.1 1.0±0.1 



TABLE VI: 0(p 6 ) fits to different sets of threshold parameters using the low energy constants h and fi. In the first row we 
fit all the threshold parameters obtained in this work. Note that the \ 2 /d.o.f. is quite large. However, as it happened in the 
equivalent fit in Sect. lIVl the largest contribution to \ comes from cp. Thus, in the following fits we include all the observables 
but cp. In the second and third rows we show the LECs obtained when cp is excluded, using / w and /o in the last term of the 
ChPT expansion, respectively. The quality of the fits notably improves. Nevertheless, the large size of the errors in the case of 
'3, fi, T2, V3 and f 4 indicates that our fits are not very sensitive to these LECs. For that reason, in the next section of the table 
( "Constrained fit to" ) we repeat the latter two fits, fixing the value of £3 to an average of lattice determinations [5( and that of 
the LECs fi to fi to the resonance saturation estimates [J. We provide an estimate of the LECs and their uncertainties as a 
weighted average of these two last fits. Let us note that the 0(p 4 ) LECs h to I4 do not lie too far from our 0(p 4 ) estimates, 
shown in the seventh row. Moreover, within our uncertainties, the resulting values of the bi parameters are very consistent 
with previous determinations, listed in the last row (I3 from fi to fi from 4] and the rest of the LECs from [ill]'). 



is also far from our 0(p 4 ) values, also given in Tabic fVTl 
or the lattice value in [5|. Similarly, those of fi to f 4 
are far from the resonance saturation estimates in 
However, all of them are still relatively compatible due 
to the resulting large uncertainties. 

As we did in the main text, we repeat this last fit to 
all parameters but cp, this time replacing / w by /o in 
the 0(p 6 ) terms, which is a change of higher order in the 
ChPT expansion. We observe that we obtain a similarly 
good description of the observables, but with LECs closer 
to the reference values from 0, [j| [ll| . 

In fact, since we observe that our fits are not very sen- 
sitive to the value of some LECs, we can ask what quality 
we can achieve if we fix these parameters to the reference 
values. We thus repeat the fit to all the observables but 
cp, fixing ^3 and f\ to to the reference values given 
in the last row of Table IVII both using /„■ and /o in the 
0(p 6 ) terms. The results of these "constrained" fits are 



shown in the fourth and fifth rows of the table. As ex- 
pected, the \ 2 /d.o.f does not increase much with respect 
to the unconstrained fit, confirming that our observables 
depend little on these LECs and thus we are not able 
to provide an estimate for them. For the rest of LECs, 
we give as estimates the weighted average of the values 
found in the constrained fits, all of them fairly compatible 
within uncertainties with the reference values. 



Appendix B: Threshold parameters in ChPT 

We show here the two-loop ChPT expressions for the 
.F-wave threshold parameters as well as the third order 
threshold parameters, ca, for all waves. The scattering 
lengths, an, and slope parameters, &£/, for the S, P and 
D waves can be found, for instance, in 
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